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Gamma-ray bursts: optical afterglows in the deep 
Newtonian phase 
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ABSTRACT 

Gamma-ray burst remnants become trans-relativistic typically in days to tens of days, 
and they enter the deep Newtonian phase in tens of days to months, during which the 
majority of shock-accelerated electrons will no longer be highly relativistic. However, 
a small portion of electrons are still accelerated to ultra-relativistic speeds and capable 
of emitting synchrotron radiation. The distribution function for electrons is re-derived 
here so that synchrotron emission from these relativistic electrons can be calculated. 
Based on the revised model, optical afterglows from both isotropic fireballs and highly 
collimated jets are studied numerically, and compared to analytical results. In the 
beamed cases, it is found that, in addition to the steepening due to the edge effect 
and the lateral expansion effect, the light curves are universally characterized by a 
flattening during the deep Newtonian phase. 

Key words: radiation mechanisms: non-thermal - stars: neutron - ISM: jets and 
outflows - gamma-rays: bursts 



1 INTRODUCTION 

Observations of X-ray, optical and radio afterglows from 
gamma-ray bursts (GRBs) definitely show that at least most 
long GRBs are of cosmological origin (e.g. Costa et al. 1997; 
Frail et al. 1997; Galama et al. 1997; Metzger et al. 1997; 
Kulkarni et al. 1998; Garcia et al. 1998; Piro et al. 1998; 
Pedersen et al. 1998; Galama et al. 1998a, b; Akerlof et al. 
1999; Vreeswijk et al. 1999; Zhu et al. 1999; Hjorth et al. 
2002). The famous fireball model, which incorporates in- 
ternal shocks to account for the main bursts, and external 
shocks to account for afterglows, becomes the most popu- 
lar model (Meszaros & Rees 1992; Rees & Meszaros 1992; 
Meszaros, Laguna & Rees 1993; Meszaros, Rees & Pap- 
athanassiou 1994; Katz 1994; Rees & Meszaros 1994; Sari, 
Narayan & Piran 1996). Observed features of GRB after- 
glows can basically be well explained in this frame-work 
(Meszaros & Rees 1997; Vietri 1997; Waxman 1997a,b; Wi- 
jers, Rees & Meszaros 1997; Meszaros, Rees & Wijers 1998; 
Sari, Piran & Narayan 1998; Dai & Lu 1998, 1999; Der- 
mer, Bottcher & Chiang 1999a; Dermer, Chiang & Bottcher 
1999b; Wijers & Galama 1999; Wijers et al. 1999). How- 
ever, we are still far from resolving the puzzle of GRBs, 
since many crucial information concerning the true nature 
of the inner engine is lost after an initial acceleration phase 
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of the fireball evolution. Studies on afterglows can hopefully 
provide important clues of GRB progenitors. For example, 
from afterglow observations, we can measure the beaming 
angle, the environment density, the total kinetic energy etc 
(Hjorth et al. 2002; Konigl & Granot 2002). For recent re- 
views on GRBs, see Piran (1999), van Paradijs, Kouveliotou 
& Wijers (2000) and Meszaros (2002). 

As the most violent bursts in the Universe since the 
Big Bang, GRBs are most impressing for their extremely 
relativistic motions, with Lorentz factor 7 > 100 — 1000. 
In 1997, Wijers et al. (1997) discussed the non-relativistic 
phase of GRB afterglows for the first time. But at the begin- 
ning of the afterglow era, the Newtonian aspects of GRBs 
were largely ignored in the literature for obvious reasons. 
The importance of non-relativistic phase was stressed and 
then extensively discussed by Huang et al. (Huang et al. 
1998a; Huang, Dai & Lu 1998b, 1999a,b; Huang et al. 2000a; 
Huang, Dai & Lu 2000b, c), who pointed out that GRB rem- 
nants enter the trans-relativistic phase (here we define it 
as 2 < 7 < 5) typically in a few to 10 days. Today, the 
importance of Newtonian phase has been realized by more 
and more authors (e.g. Kobayashi, Piran & Sari 1999; Livio 
& Waxman 2000; Chevalier & Li 2000; Frail, Waxman & 
Kulkarni 2000; Weiler et al. 2002; Kulkarni 2002). 

Afterglows from both isotropic fireballs and highly col- 
limated jets in the trans-relativistic and non-relativistic 
phases have been investigated in great detail by Huang et 
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al. (1998b, 1999a,b, 2000a,b,c). Although their dynamical 
equations are valid for sufficiently long period, their dis- 
cussion cannot extend to the deep Newtonian phase, dur- 
ing which most shock-accelerated electrons will cease to be 
ultra-relativistic and will no longer emit synchrotron radia- 
tion. In this article, a more reasonable distribution function 
will be derived for electrons, so that we can extend our cal- 
culation into the deep Newtonian phase. The structure of 
our article is arranged as follows. We describe our dynami- 
cal and radiation model in Section 2. Based on the revised 
model, optical afterglows from isotropic fireballs, conical jets 
as well as cylindrical jets are then investigated numerically 
in Section 3. Emphases will be put on the light curve be- 
haviour in the deep Newtonian phase. Section 4 is a brief 
discussion. 



2 MODEL 

The overall dynamical evolution of GRB remnants has been 
studied by many authors (Huang et al. 1999a,b; Kobayashi 
ct al. 1999; Dermer & Humi 2001; Panaitescu & Kumar 
2001a, 2002). Here we adopt the simple dynamical model 
proposed by Huang et al. (1999a,b), which is characterized 
mainly by the following differential equation, 



d 7 
dm 



7 2 "1 



M ej + em + 2(1 - e)^m '' 



(1) 



where 7 is the bulk Lorentz factor of the shocked interstellar 
medium (ISM), M c j is the initial ejecta mass, e is the radia- 
tive efficiency and m is the swept-up ISM mass. It has been 
shown that this equation is correct for both ultra-relativistic 
and non-relativistic blastwaves, no matter whether they are 
adiabatic or highly radiative (Huang et al. 1999a,b). In re- 
alistic GRB remnants, the blastwaves can be radiative only 
within the initial two or three hours following the main burst 
(Dai, Huang & Lu 1999). They will become completely adi- 
abatic after that period. So, for simplicity, we will take e = 
directly in this study. 

For beamed ejecta with a half opening angle of 6, Huang 
et al. (2000a, b) have proposed an improved method to de- 
scribe the lateral expansion of the remnant. Their descrip- 
tion is based on a refined expression for the sound speed 
that is appropriate in both relativistic and non-relativistic 
phases. In Section 3, we will solve the dynamical evolution 
of isotropic fireballs, conical jets as well as cylindrical jets 
numerically. For details of the numerical procedure, we refer 
readers to Huang et al. (1999b, 2000a,b) and Cheng, Huang 
& Lu (2001). 

Synchrotron radiation from the shock-accelerated ISM 
electrons plays a major role in the optical afterglows. To 
calculate this emission, two factors are essential to be known: 
the magnetic field and the energy of electrons. As usual, we 
assume that the magnetic energy density in the comoving 
frame is a fraction £| of the total thermal energy density 
(B' 2 /8ir — £|e'). Electrons are also usually assumed to carry 
a fraction £ e of the proton energy and follow a power-law 
distribution according to their Lorentz factors, 



dK 

d-) c 



OC 7 C P , (7e,min < 7e < 7e,max), 



(2) 



m c p — 1 

with m p and m c being the proton and electron mass re- 
spectively. The maximum Lorentz factor is constrained by 
synchrotron cooling, 7 c , m ax ~ 10 8 (-B'/1G) _1/ ^ 2 (Meszaros et 
al. 1993; Vietri 1997; Totani 1999). Sari et al. (1998) fur- 
ther suggested that radiation loss may play an important 
role in the process and can change the distribution function 



-(p+i) 



for electrons above a critical Lorentz 



to dNc/d'jc oc 7< 
factor 7c, which is given by j c = 6Tvm c c/(aT'fB' 2 t), where 
or is the Thomson cross section and t is the time of the 
observer. Based on these considerations, a detailed descrip- 
tion of electron distribution has been presented by Dai et al. 
(1999, also see Huang et al. 2000a,b). 

The above discussion is valid only when electrons are 
ultra-relativistic, i.e., 7 c ,mm 3> 1. However, in the deep New- 
tonian phase, this condition will no longer be satisfied. The 
bulk Lorentz factor (7) will be so close to unity that 7 c ,mm 
will also be very close to one. It means the majority of elec- 
trons will cease to be ultra-relativistic. For example, taking 
typical values of £ e = 0.1 and p = 2.2, 7 c ,mm will equal to 
5 when 7 fa 1.13 (or (3 — ^f 1 — 1/7 w 0.47). Anyway, the 
distribution of electrons with respect to their kinetic ener- 
gies can still be assumed to be a power-law function, which 
now takes the following form, 



OC (70 - 1) P , (7e ,min < 7c < 7e,max). 



(4) 



dm 

d7e 

Most electrons are now non-relativistic and their cyclotron 
radiation cannot be observed in optical bands. But there are 
still many relativistic electrons with Lorentz factor above 
a critical value, 7o, sy n, and still capable of emitting syn- 
chrotron radiation. In this case, to calculate the optical af- 
terglow, we just need to integrate the emission over all those 
electrons satisfying 7c,syn < 7c < 7c, max- The determination 
of 7c,syn is somewhat arbitrary. We believe that it could be 
some value like 5 or 10. Anyway, it is lucky that this un- 
certainty will not bring any essential effect into the optical 
afterglows (see Section 3). 

Based on equation (4) , we can now re-derive the electron 
distribution as follows, following Dai et al.'s (1999) treat- 
ment (also see Huang et al. 2000a, b) , 

(i) For 7 C < 7c 

,min 5 

^ =d(7c-ir (P+1) , (7c,min < 7c < 7c,max) , (5) 
07e 



(7c, mill - 1) P - (7c, max - 1)~ 



-N ei 



(6) 



where A^ic is the total number of radiating electrons in- 
volved. 

(ii) For 7c ,min < 7c < 7c ,raax 7 



dK f C 2 (7c-l)- p : 



d 7 c \ Cs(7e-1) 
where 

C 2 = C 3 /( 7 c - 1) , 



-(p+i) 



(7c ,min < 7c < 7c), 

(7c < 7c < 7c, max), 



(7) 



(8) 



C, 



where the minimum Lorentz factor can be easily derived as, 



(7e 



l) 1 "* - (7c - I) 1 " 



(p-l)( 7c -l) 



© 2003 RAS, MNRAS 000, 000-000 



3 



(7c -l)- p - (7c, max - l)- p 



P 



N, 



clc • 



(iii) If 7c > 7c,max, then 

(7e 

where 

C 4 



dN' 



< 7c < 7c 



p-1 



(7c,min - l) 1 "? - (7e,max - l) 1 "" 



N ele . 



(9) 



(10) 



(11) 



Of course, 7c 3> 1 and 7 c ,max 3> 1 are safely satisfied, so 
that 7c — 1 and 7 c ,max — 1 can be simplified as j c and 7 c ,max 
respectively in these equations. 



been detected for more than six months, and radio after- 
glows are even detectable one or three years later (Frail et 
al. 2000; Kulkarni 2002), we see that the study of afterglows 
in the deep Newtonian phase is really essential. 

In our calculations, we have assumed that the mini- 
mum Lorentz factor of electrons capable of emitting syn- 
chrotron radiation is 7 c , S yn = 5. This evaluation is some- 
what arbitrary, but it does not affect the optical light curve 
too much. The characteristic synchrotron frequency of an 
electron with Lorentz factor 7 C in a magnetic field B' is 
v = 7 c 2 eB'/(27rm c c) w 2.8 x 10 6 7c (B'/1G) Hz. Electrons in 
the lower energy range thus do not contribute to optical flux 
density significantly. In fact, we have taken 7e,syn as 10 or 
even 50 in our trial calculations and just found no difference 
in the optical light curves. 



3 NUMERICAL RESULTS 

Afterglows from isotropic fireballs as well as collimated jets 
have been studied in great detail by Huang et al. (Huang 
ct al. 1999a,b, 2000a,b,c; Cheng, Huang & Lu 2001). How- 
ever, those discussion in general can not extend to the deep 
Newtonian phase. In this section, we use our revised model 
to continue the study. In all our calculations here, we take 
the following typical parameter values: 7 c ,syn = 5, £ c = 
0.1, .Dl = 1 Gpc and p = 2.2, where Dl is the luminosity 
distance. The initial value of 7 is fixedly taken as 70 = 300. 
Other parameters, such as the total isotropic kinetic energy 
of the blastwave (Eo), the number density of the interstellar 
medium (n) and the magnetic energy fraction (£|), will be 
given separately elsewhere. The effects of equal arrival time 
surfaces (Waxman 1997c; Sari 1998; Panaitescu & Meszaros 
1998; Huang et al. 2000a,b) are taken into account in our 
studies. For more details of the calculation, readers are ref- 
ered to Huang et al. (1999b, 2000a,b) and Cheng, Huang & 
Lu (2001). 



3.1 Isotropic fireballs 

Fig. 1 shows the evolution of the Lorentz factor for some 
exemplary isotropic fireballs and the corresponding _R-band 
flux density (Sr). Analytical solution requires that the fire- 
ball ev olves a s 7 oc i~ 3//s in the ultra-relativistic phase and 
P — yj"{ 2 — I/7 cc t~~ 3 / 5 in the Newtonian limit (Wijers et 
al. 1997; Huang et al. 1999a, b). Fig. la consists with the re- 
quirement exactly. Wijers et al. (1997) and Dai & Lu (1999, 
2000) have also derived the theoretical afterglow light curve 
analytically for both relativistic and non-relativistic phases. 
In case of slow cooling, we have 



Sr oc 



t (3-3p)/4 
£(21-15p)/10 



(7 » 1), 

03 « 1). 



(12) 



Taking p = 2.2, we get Sr oc t~ - 9 and Sr oc i -1 ' 2 for rel- 
ativistic and Newtonian phases respectively. Our numerical 
results are consistent with their solutions. For example, the 
slope of the solid line in Fig. lb is ~ —0.92 and ~ —1.22 in 
the relativistic phase and the non-relativistic phase respec- 
tively. 

Fig. 1 also clearly shows that the remnant enters the 
deep Newtonian phase in a relatively short period. For ex- 
ample, 7c,min becomes less than 7 c , S yn at ~ 10 7 s, i.e., about 
3 months. Since optical afterglows from some GRBs have 



3.2 Conical jets 

Collimation is very important in GRBs, which can affect 
the intrinsic kinetic energy and may provide direct clues 
on the progenitors (Frail et al. 2001; Ioka & Nakamura 
2001; Panaitescu & Kumar 2001a,b, 2002; Ramirez-Ruiz & 
Lloyd- Ronning 2002; Yamazaki, Ioka & Nakamura 2002). 
Beaming effects in afterglows have been discussed exten- 
sively in the literature. It is generally believed that due to 
both the edge effect (Kulkarni et al. 1999; Meszaros & Rees 
1999; Panaitescu & Meszaros 1999) and the lateral expan- 
sion effect (Rhoads 1997, 1999), afterglows from a conical jet 
are characterized by a break in the light curve. The break 
point is approximately determined by 7 ~ I/60, where 80 
is the initial half opening angle. GRBs 990123, 990510 and 
000301c are regarded as good examples (Castro-Tirado et 
al. 1999; Harrison et al. 1999; Kulkarni et al. 1999; Sari, 
Piran & Halpern 1999; Wijers et al. 1999; Masetti et al. 
2000; Holland et al. 2000; Berger et al. 2000; Jensen et al. 
2001; Huang et al. 2000c). However, detailed numerical stud- 
ies show that the break is usually quite smooth (Moderski, 
Sikora & Bulik 2000). Huang et al. further found that the 
break is parameter-dependent (Huang et al. 2000a, b) . They 
also suggested that the trans-relativistic and non-relativistic 
stages are important and should be considered carefully. In 
this section, we use our revised model to reveal the behaviour 
of conical jets in the deep Newtonian phase. For more details 
of the calculation, please refer to Huang et al. (2000a,b). 



3. 2. 1 Conical jets without lateral expansion 

Fig. 2 illustrates the evolution and afterglows of some exem- 
plary conical jets without lateral expansion. In fact, the dy- 
namical evolution of such a jet should be very similar to that 
of an isotropic fireball, i.e., 7 oc t~ :i ' 8 in the ultra-relativistic 
phase and (3 oc t~ 3/l5 in the Newtonian stage. Fig. 2a shows 
these trends exactly. As for the optical light curve, analyt- 
ical solutions predict that it should follow Sr oc ^ 3_3p )/ 4 
(t~°- 9 for p = 2.2) before 7 = l/0 = 10. But soon after 
the 7 = 10 point, a break should appear and the light curve 
becomes Sr oc i~ 3p/4 (t" 1 ' 65 for p = 2.2). This relation can 
be derived as follows. Analytically we have 7 oc i~ 3//g and 
R oc t 1 / 4 . Then the characteristic synchrotron frequency is 
i/ m oc 77c,mi n S' oc 7 4 oc t~ 3 ^ 2 , and the corresponding maxi- 
mum flux density is S„, miK oc N e i e -yB' /-y~ 2 oc R 3, y 4 oc t~ 3 / 4 , 
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where the 7 2 item in S„, mM arises from the edge effect 
when 7 < 1/8q. So, the observed optical flux density is 



Sr « S lAinax 



"(p-1)/2 



OC t 



-3p/4 



Such a light curve 



break can be clearly seen in Fig. 2b. For example, the solid 
line can be fit as Sr oc t~ - 89 during 10 3 — 10 4 s and it is 
approximately Sr. oc t~ 180 during 10 5 — fO 6 s. The slopes 
are acceptably consistent with analytic values. However, it 
should be noted that the break occurs obviously later than 
the time determined by 7 = l/9 — 10, which has been 
pointed out by Huang et al. in previous studies (Huang et 
al. 2000a,b). 

Maybe the most striking feature in Fig. 2b is the flatten- 
ing of the light curves in the deep Newtonian phase. This is 
not difficult to understand. In such a non-relativistic phase, 
the remnant, although beamed, should behaves more or less 
like a spherical shell since the edge effect no longer exists 
(Livio & Waxman 2000) . Beaming will only reduce the total 
flux density by a constant factor. In fact, the slope of the 
solid line is approximately —1.2 after 10 9 s, which is just 
the value expected for an isotropic fireball (see Section 3.1). 
So, we see that the overall evolution of optical afterglows 
from a conical jet without lateral expansion can be typically 
described by, 



j(3-3p)/4 ^ 
£(21-15p)/10 



(7 > l/0o), 

(7 < l/0o and/3 

09 « 1)- 



(13) 



3.2.2 Conical jets with lateral expansion 

In realistic cases, jets may expand laterally at comoving 
sound speed. Fig. 3a shows the evolution of Lorentz factor 
for such realistic conical jets. The curves can be approxi- 
mately fit by 7 oc t -0 ' 39 in the highly relativistic phase and 



by 7 — 1 oc t 



in the non-relativistic phase, consistent 



1.03 
0.14 



with the analytic results of 7 oc t 3 ^ 4 and 7 — 1 oc t 6 ^ 5 
respectively. Again we see that 7 e ,min becomes less than 5 
when t > 3 x 10 6 — 2 x 10 7 s. 

Fig. 3b illustrates the i?-band afterglows correspond- 
ingly. We see that the general behaviour of the light curves 
is very similar to that in Fig. 2b. However, slight difference 
can still be found and deserves addressing in some detail. For 
example, the solid line approximately follows Sr. oc t 
in the ultra-relativistic phase, which is steeper by ~ t 
than that in Fig. 2b. It indicates that the lateral expansion 
tends to make the afterglow decay faster. Subsequently the 
light curve steepens to Sr oc i~ 2 ' 24 , which is also steeper by 
~ i 0,44 than that in Fig. 2b. It is clear that lateral expan- 
sion contributes to the light curve break significantly. In the 
deep Newtonian phase, the light curves are again markedly 
characterized by a flattening. The slope of the solid line in 
this segment is ~ —1.26, also consistent with the analytical 
result of —1.2 for an isotropic fireball. In short, we conclude 
that the overall evolution of optical afterglows from a conical 
jet with lateral expansion can be typically expressed as, 



Sr oc 



t (3-3p)/4 
£(21-15p)/10 



(7 > 1/flo), 

(7< l/0o and/3~ 1), 

09 « 1)- 



(14) 



3.3 Cylindrical jets 

The geometry of GRB jets is usually assumed to be conical. 
However, many of the relativistic outflows in radio galaxies 
are found to maintain constant cross-sections at large scales 
(Perley, Bridle & Willis 1984; Biretta, Sparks & Macchetto 
1999). This has led Cheng et al. (2001) to suggest that GRB 
jets might also be cylindrical. In fact, a very similar idea, i.e., 
the so called cannon-ball model, has been suggested as the 
GRB trigger mechanism by Dar et al. (Shaviv & Dar 1995; 
Dar 1998), and observed afterglows have also been examined 
carefully in this frame (Dado, Dar & De Rujula 2002). 

Theoretical GRB afterglows from cylindrical jets have 
been studied by Cheng et al. (2001) in great detail. However, 
their discussion again cannot extend to the deep Newtonian 
phase. Here we use our revised model to continue the study. 
Since a cylindrical jet without lateral expansion generally 
decelerates very slowly and can remain to be highly rela- 
tivistic for as long as ~ 10 8 — 10 9 s (Cheng et al. 2001), we 
will discuss only cylindrical jets with lateral expansion. 

Fig. 4a shows the evolution of 7 for some exemplary 
cylindrical jets. In the highly relativistic phase the curves 
approximately follow 7 oc t~ ' 5 , and in the non-relativistic 
phase they follow 7 — 1 oc i~ 1,17 . The timing indices are 
consistent with the analytical results of —1/2 and — 6/5 
respectively (Cheng et al. 2001). Note that in this figure, 
7c,min is already less than 5 after ~ (2 — 4) x 10 6 s. Fig. 
4b illustrates the evolution of optical afterglows correspond- 
ingly. In the ultra-relativistic phase, the solid line can be 
approximated as Sr oc i~ 2,34 , consisting with the analyti- 
cal solution of Sr oc t~ p (Cheng et al. 2001). After entering 
the non-relativistic phase, the remnant is expected to re- 
semble an isotropic fireball, so that the light curve should 
be Sr oc t (21 - 15p)/1 ° (Cheng et al. 2001). Taking p = 2.2, 
we get the analytical timing index of —1.2. In fact, such a 
flattening is really observed in Fig. 4b, and the slope of the 
solid line is just ~ —1.33 in the deep Newtonian phase. From 
these studies, we see that the optical afterglow from laterally 
expanding cylindrical jets can typically be well represented 
by, 



Sr OC 



(21-15p)/10 



(7 » 1), 

(/3 « 1). 



(15) 



4 DISCUSSION AND CONCLUSIONS 

We have shown that typically in ~ 10 7 s, the majority of 
shock-accelerated electrons will become non-relativistic, i.e., 
7o,min < 7o,syn- These electrons can only emit cyclotron ra- 
diation which falls mainly in the frequency range of ~ 10 
- 10 7 Hz. However, a small portion of electrons are still 
highly relativistic and capable of emitting synchrotron radi- 
ation. A revised distribution function has been derived for 
electrons, which allows us to calculate the emission from 
these relativistic particles. 

Optical afterglows from both isotropic fireballs and 
highly collimated jets are studied numerically based on the 
revised model. For isotropic fireballs, the optical light curve 
steepens by i( 15 p~ 27 )/ 20 after entering the deep Newtonian 
phase. This corresponds to ~ t for p — 2.2, or ~ t for 
p = 3. For conical jets, we find that the light curve steepens 
slightly later than the moment determined by 7 = l/0o- It 
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is clear that this steepening is due to the edge effect and 
is strengthened by lateral expansion. Additionally, the op- 
tical light curves are universally characterized by a striking 
flattening during the deep Newtonian phase, which can be 
attributed to the resembling of the remnant to an isotropic 
fireball at such late stages (Livio & Waxman 2000). The 
flattening is also present in the light curves of cylindrical 
jets. 

Long-lasting optical afterglows have been observed from 
a number of GRBs, e.g. GRBs 970228, 970508, 991208 
(Fruchter et al. 1999, 2000; Castro- Tirado et al. 2001). Op- 
tical transients in these events are typically detectable for 
~ 100 — 200 d. Our studies here are of obvious importance 
in such cases. However, the flattening of the light curve in 
the deep Newtonian phase predicted for highly collimated 
jets has not been observed in realistic observations of jet 
candidates. A possible reason may be that the optical flux 
densities are usually too faint. In fact, to observe this flat- 
tening clearly, we should follow the optical transient for at 
least ~ 10 8 s, i.e., more than ~ 3 years. This could hardly 
be actualized even with the Keck telescope and the Hub- 
ble Space Telescope. Additionally, optical emission from the 
host galaxy makes the measurement even more difficult. 
Anyway, the next generation optical telescopes may cast 
some light on the problem. 

Radio afterglows can hopefully be observed for even 
longer period. In fact, a 450-day light curve of the radio 
afterglow of GRB 970508 has been reported by Frail et al. 
(2000). GRB 980703 has even been monitored in radio bands 
for more than 1000 days (Kulkarni 2002). At such stages, 
the GRB remnant should have entered the deep Newtonian 
phase. Observations of radio afterglows on such a late stage 
have the priority that they allow for a direct measure of 
the intrinsic kinetic energy (Frail et al. 2000). The effect of 
7c,min will have to be considered in fitting those radio light 
curves theoretically. Calculation of radio emission involves 
synchrotron self absorption, which is beyond the scope of 
this article and will be discussed elsewhere. 

Orphan afterglow survey is thought to be a useful 
method to determine the degree of beaming in GRBs 
(Rhoads 1997; Meszaros, Rees & Wijers 1999; Totani & 
Panaitescu 2002; Nakar, Piran & Granot 2002; Granot et 
al. 2002). However, the method is also troubled by many 
problems (Dalai, Griest & Pruet 2002; Huang, Dai & Lu 
2002; Levinson et al. 2002; Nakar & Piran 2003). For ex- 
ample, as suggested by Huang et al. (2002), the method is 
greatly complicated by the possibility that there might be 
many failed GRBs, i.e., isotropic fireballs with 1 <C 7 <C 100, 
which cannot successfully generate 7-ray bursts but do can 
produce afterglows. Anyway, a successful application of the 
method would involve two important ingredients: the orphan 
needs to be discovered as early as possible and needs to be 
followed as long as possible. In fact, it has been suggested 
that radio searches, which are capable of finding potential 
GRB remnants as old as hundreds of days, may provide the 
best hope to find the missing orphans (Dalai et al. 2002; 
Totani & Panaitescu 2002; Levinson et al. 2002). Again we 
see that theoretical investigation of GRB afterglows in the 
deep Newtonian phase is essential. 
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Figure 1. Evolution of the Lorentz factor for isotropic fireballs 
(a) and the corresponding i?-band optical afterglows (b) . We take 
£j| = 10~ 4 . Eq and n values are marked in units of erg and 
cm -3 respectively. Other parameters are evaluated in the first 
paragraph of Section 3. The black dot on each curve indicates 
the moment when 7 Cjm i n = 7o,syn, and the open circles on the 
light curves mark the time when 7 = 2,5 and 10 respectively. 
Note that the evolution of 7 is almost identical for the cases of 
Eq = 10 52 erg, n = 10 cm~ 3 and Eq = 10 51 erg, n = 1 cm~ 3 . 
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Figure 2. Evolution of the Lorentz factor for conical jets without 
lateral expansion (a) and the corresponding optical afterglows 
(b). We take 9q = 0.1, §j| = 10 -4 . Isotropic equivalent energy Eq 
and n values are marked in units of erg and cm -3 respectively. 
Other parameters are evaluated in the first paragraph of Section 
3. Observers arc assumed to be on the axis of the jet. The black 
dots and the open circles have the same meaning as in Fig. 1. 
Note that the evolution of 7 is almost identical for the cases of 
Eq = 10 53 erg, n = 10 cm~ 3 and Eq = 10 52 erg, n = 1 cm~ 3 . 



Figure 3. Same as Fig. 2 except that the conical jets here are 
expanding laterally at comoving sound speed. 
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Figure 4. Evolution of the Lorentz factor for cylindrical jets with 
lateral expansion (a) and the corresponding optical afterglows (b). 
We take £| = 0.01 and assume that the initial comoving lateral 
radius of the jet is ao = 0.01-Rq. Isotropic equivalent energy Eg 
and number density n are marked in units of erg and cm -3 re- 
spectively. Other parameters are evaluated in the first paragraph 
of Section 3. Observers are assumed to be on the axis of the jet. 
The black dots and the open circles have the same meaning as 
in Fig. 1. Note that the evolution of 7 is almost identical for the 
cases of E = 10 54 erg, n = 100 cm -3 and E = 10 53 erg, n = 10 
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